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We determine the strong-isospin violating component of the neutron-proton mass difference from 
fully-dynamical lattice QCD and partially-quenched QCD calculations of the nucleon mass, con- 
strained by partially-quenched chiral perturbation theory at one-loop level. The lattice calcula- 
tions were performed with domain-wall valence quarks on MILC lattices with rooted staggered 
sea-quarks at a lattice spacing of b ~ 0. 125 fm, lattice spatial size of L ~ 2.5 fm and pion masses 
ranging from m K ~ 290 MeV to ~ 350 MeV. At the physical value of the pion mass, we predict 
M„ - M p | ^ = 2.26 ± 0.57 ± 0.42 ± 0. 1 MeV where the first error is statistical, the second error 
is due to the uncertainty in the ratio of light-quark masses, rj = m„/wrf, determined by MILC [[[J, 
and the third error is an estimate of the systematic due to chiral extrapolation. 
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1. Introduction 

It is a basic property of our universe that the neutron is slightly more massive than the pro- 
ton. The electroweak interactions are responsible for this mass difference, which receives con- 
tributions from two sources. The strong-isospin breaking contribution (also known as charge- 
symmetry breaking, for a review see Ref. |j2p) is due to the difference in the masses of the up 
and down quarks, ultimately determined by the values of the Yukawa couplings in the Standard 
Model of electroweak interactions and the vacuum expectation value of the Higgs field. With 
the down-quark more massive than the up-quark, the neutron would be more massive than the 
proton in the absence of electromagnetism. The experimental neutron-proton mass difference of 
M n -M p = 1.2933317 ±0.0000005 MeV ^ receives an estimated electromagnetic contribution 
of 0] M n -M p \ em = -0.76 ±0.30 MeV (estimated using the Cottingham sum-rule saturated by 
the Born diagrams), and the remaining mass difference of M n - M p | " = 2.05 T 0.30 MeV is due 
to strong-isospin breaking. The uncertainty is estimated from the Born contribution of baryon- 
resonances to the Cottingham sum-rule. 

Before delving into the details of the lattice calculation we have performed, it is useful to 
remind ourselves why we should care about this mass-splitting in the first place. If we were in a 
situation where M n < M p +M e — M v , hydrogen would not be stable. It would decay weakly via 
p + e — > n ± v and chemistry of our universe would be quite different from what we are familiar 
with. Therefore, in the limit of strong-isospin symmetry, hydrogen would not be stable, and from 
our present calculation we find that for hydrogen to be stable, m„/mrf< 0.77 (taken from the central 
values of our determination and the electromagnetic contribution). Of course, one would still have 
some form of nuclear physics, but it might look quite different from what we are familiar with. In 
order to determine just how similar, or different, nuclear physics would be, we need to perform 
further calculations. 



2. The chiral expansion of the Nucleon Mass 



In order to extract isospin breaking quantities from lattice calculations on isospin symmetric 
lattices, one must perform partially-quenched calculations, and then determine coefficients in the 
chiral theory describing QCD quantities. In QCD the proton mass has a chiral expansion of the 
form 
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where tj = m M /m^ is the ratio of the mass of the up-quark and down-quark. The function F K is 
given by F n = F{m n , A,/i), with 
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and arises from A-resonance intermediate states in one-loop diagrams. A is the mass difference 
between the A-resonance and the nucleon. The mass-splitting between the neutron and proton due 
to the quark mass differences is therefore, 

M n -M p \ d - = ^(2«-£)(^) "4 • (2-3) 

The one-loop contributions at ff(m 3 J 2 ) cancel in the mass-difference, as the pions are degenerate 
up to 0(m 2 q ). Consequently, at leading order, it is the combination of parameters 2a — j8 that needs 
to be extracted from lattice calculations. 

The nucleon mass was computed in SU (2)l <8> SU (2)r partially-quenched chiral perturbation 
theory (PQ^PT) in Ref. [Q], and involves quite lengthy expressions. In the case of an isospin- 
symmetric set of lattice configurations, the proton masses for different combinations of valence 
quarks can be found in Ref. As an example, the mass of a proton with valence quantum 
numbers Vi , Vi , V2 on a sea of V\ quarks has the form 
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where the function S is given by S n = S(m n ,A,n) with 
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m v r 8 ,Vi;V r c is the mass of a pion composed of valence quarks V a , Vb on configurations with quarks Ve- 
in partially-quenched proton masses, there are contributions that depend upon the axial-isosinglet 
coupling, g\. This coupling is one of the fit parameters in our extraction. 

3. The Lattice QCD Calculation 

Our computation [|]] uses the mixed-action lattice QCD scheme developed by LHPC [^, using 
domain-wall valence quarks from a smeared-source on Nf = 2 + 1 asqtad-improved [Q, [T(J MILC 



configurations generated with rooted staggered sea quarks [ |18fl that are hypercubic-smeared 



'For recent discussions of the "legality" of the mixed-action and rooting procedures, see Ref. JlH |l2[ |l3[ |l4[ |l5[ 
[l^ , [Tt[ ] . While it still could be the case that staggered fermions do not reproduce QCD in the continuum limit, recent 
work by Shamir Jl7| ] on this issue, relying heavily on the renormalization group evolution of the staggered theory, has 
concluded that under certain plausible assumptions that can be tested, the continuum limit of the staggered theory is 
QCD. 
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Ensemble 


Theory 


bnii 


bm dwf 


10 3 x bm res 


# props. 


2064f21b679m007m050 


QCD 


0.007 (Vj) 


0.0081 (Vi) 


1.604 ±0.038 


468x3 


2064f21b679m007m050 


PQ 


0.007 (Vi) 


0.0138 (V 2 ) 


1.604 ±0.038 


367x3 


2064f21b679m007m050 


PQ 


0.007 (Fi) 


0.0100 (V 3 ) 


1.604 ±0.038 


367x2 


2064f21b679m010m050 


QCD 


0.010 (V 2 ) 


0.0138 (V 2 ) 


1.552 ±0.027 


658x3 


2064f21b679m010m050 


PQ 


0.010 (V 2 ) 


0.0081 (Vj) 


1.552 ±0.027 


658x1 



Table 1: The parameters of the MILC gauge configurations and domain-wall propagators used in this work. 
For each propagator the extent of the fifth dimension is L5 = 16. The notation of quarks, Vi,V2,V3, is defined 
in the text. The last column is the number of propagators generated, and corresponds to the number of lattices 
times the number of different locations of sources on each lattice. We have worked at one lattice-spacing 
only, b ~ 0. 125 fm. PQ denotes partially-quenched QCD. For each set of configurations bm s — 0.050. 



(HYP-smeared) [ |19| , |20[ |21|, |22j]. In the generation of the MILC configurations, the strange-quark 
mass was fixed near its physical value, bm s = 0.050, (where b ~ 0. 125 fm is the lattice spacing) 
determined by the mass of hadrons containing strange quarks. The two light quarks in the con- 
figurations are degenerate (isospin-symmetric). The domain-wall height is m = 1.7 and the extent 
of the extra dimension is L5 = 16. The MILC lattices were "chopped" using a Dirichlet bound- 
ary condition from 64 to 32 time-slices to save time in propagator generation. In order to extract 
the terms in the mass expansion, we computed a number of sets of propagators corresponding to 
different valence quark masses, as shown in Table I. 

On 468 bm; = 0.007 (denoted by V\) lattices we have computed three sets corresponding to 
the QCD point with a valence-quark mass of bmdwf = 0.0081 (Vi), three sets on 367 bm\ = 0.007 
lattices with a valence quark mass of bmdwf = 0.0138 (denoted by V%), and two sets with a valence 
quark mass bmdwf = 0.0100 (denoted by V3). On 658 of the bm\ = 0.010 (V 2 ) lattices we have 
computed three sets at the QCD point with a valence-quark mass of bmdwf = 0.0138 (V 2 ) and 
one set with a valence quark mass of bm dw f = 0.0081 (V\). The parameters used to generate 
the QCD-point light-quark propagators have been "matched" to those used to generate the MILC 
configurations so that the mass of the pion computed with the domain-wall propagators is equal (to 
few-percent precision) to that of the lightest staggered pion computed with the same parameters as 



the gauge configurations [18]. The lattice calculations were performed with the Chroma software 
suite [^3|, 24] on the high-performance computing systems at the Jefferson Laboratory (JLab). To 
extract the isospin-breaking terms it is sufficient to look at mass-differences between proton states. 
The ratio of correlation functions is formed so that a mass-difference can be extracted from its 
large-time behavior. An example of the ratio of correlation functions found in our work is shown 
in fig. [j] We define the mass difference to be 

AM p (V a ,V b ,V c ;V d ) = M p (V a ,V b ,V c ;V d ) - M p (V d ,V d ,V d ;V d ) . (3.1) 

Each pair of partially-quenched propagators generated four different proton states and hence 
three different mass-splittings from the isospin symmetric QCD proton state. The results of our 
calculations are shown in fig. |[ By fitting four parameters, a,f},gi and gAN to the six mass- 
splitting we are able to extract the combination of constants in eq. ^3|. Further, using the ratio of 
light-quark masses determined by the MILC collaboration, m u /m d = 0.43 ±0.01 ±0.08 [[[]], we can 
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Figure 1: The effective mass plot for the mass splitting AM p (V\,Vi_,V2',V\). 
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Figure 2: 77;e partially-quenched proton mass differences (in MeV) calculated from the bmi — 0.007 and 
0.010 MILC lattices plotted vs the pion mass composed of sea quarks. Various data have been displaced 
horizontally by small amounts for display purposes. A lattice spacing ofb = 0.125 fm has been used to set 
the scale. 



make a prediction for the strong-isospin breaking contribution to the neutron-proton mass-splitting 
of 



M„-M, 



\d—u 



2.26 ±0.57 ±0.42 ±0.10 MeV , 



(3.2) 



where the last error is an estimate of the systematic error due to truncation of the chiral expansion. 
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This is to be compared with a value of M n — M p \ " = 2.05 =F 0-30 MeV as determined via the 
Cottingham sum-rule. 

4. Conclusions 

We have performed partially-quenched calculations of the proton mass with domain-wall va- 
lence quarks on the isospin-symmetric coarse staggered MILC lattices [H). These calculations have 
allowed us to isolate the contribution to the neutron-proton mass-splitting due to the light-quark 
mass-splitting at leading and next-to-leading order in the SU (2)l ® SU (2)r chiral expansion. We 
find a value that is consistent with the estimate arrived at from the experimental mass-splitting and 
the best estimate of the electromagnetic conuibution. It is clear that with more computer power, 
this important quantity will be calculated to high precision. 
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